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We present the first fully ab initio calculations for attosecond transient absorption spectroscopy of 
atomic krypton with overlapping pump and probe pulses. Within the time-dependent configuration 
interaction singles (TDCIS) approach, we describe the pump step (strong-field ionization using 
a near-infrared pulse) as well as the probe step (resonant electron excitation using an extreme- 
ultraviolet pulse) from first principles. We extent our TDCIS model and account for the spin-orbit 
splitting of the occupied orbitals. We discuss the spectral features seen in a recent attosecond 
transient absorption experiment [A. Wirth et al., Science 334, 195 (2011)]. Our results support 
the concept that the transient absorption signal can be directly related to the instantaneous hole 
population even during the ionizing pump pulse. Furthermore, we find strong deformations in the 
absorption lines when the overlap of pump and probe pulses is maximum. These deformations can 
be described by relative phase shifts in the oscillating ionic dipole. We discuss possible mechanisms 
contributing to these phase shifts. Our finding suggests that the non-perturbative laser dressing of 
the entire iV-electron wave function is the main contributor. 

PACS numbers: 32.80.Rm,42.65.Re,31.15.A- 
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I. INTRODUCTION 

The interaction of matter with light is a key process in 
physical systems on any length scale. The fundamentals 
of matter-light interaction can be best studied in atomic 
systems due to their relative simplicity. The absorption 
of light promotes electrons into excited states. If enough 
energy is absorbed by the system, one or more electrons 
can leave the atom, i.e., ionization takes place [1-3]. 
The most common types of ionization are: single-photon 
and few-photon ionizations [4-6] , above-threshold ioniza- 
tion [7-11], and tunnel ionization [12-17]. 

Recently, high-harmonic generation (HHG) has be- 
come a major tool in attosecond physics, allowing one to 
generate ultrashort light pulses with broad spectral band- 
widths [18, 19]. From the ability to generate attosecond 
pulses [20], an entire new research area has emerged [21] 
focusing on electronic dynamics [22-28] and molecular 
motion [29-32] on their fundamental time scale. A par- 
ticularly interesting aspect is the electron motion and 
the corresponding hole creation dynamics during the ion- 
ization process [33-35]. The high pulse intensities used 
in these experiments distort significantly the potential 
of the electrons such that it is possible for the electron 
to tunnel through or even travel over the barrier out of 
the system (i.e., tunnel ionization or barrier suppression 
regime, respectively). 

A well-known model to describe tunnel ionization in 
atomic systems is the Ammosov, Delone, and Krainov 
(ADK) model [36, 37]. It applies to intense low-frequency 
fields, where a quasistatic approximation can be made, 
meaning electrons follow adiabatically the external field. 
In these kind of fields, tunneling rates and final ion pop- 
ulations can be well reproduced by the ADK model [38]. 



Short, few-cycle pulses give access to instantaneous 
rather than cycle-averaged quantities and reveal the non- 
adiabatic behavior of the electronic motion. In this case, 
an explicit time-dependent treatment of the ionization 
process is advantageous. In combination with a multi- 
channel theory, the dynamics in the relative phases be- 
tween generated ionic states can be captured; something 
that cannot be done by the ADK model. In the past, the 
state of the ion (after ionization) has not been of high 
interest, since it was not experimentally accessible. 

Such a technique does now exist: attosecond transient 
absorption spectroscopy [39-42]. Transient absorption 
spectroscopy has been used for years to study chemical 
reactions on the femtosecond time scale [43]. However, 
just recently this technique has been extended to the at- 
tosecond regime [39], where it is possible to probe the 
diagonal and off-diagonal elements of the ion density ma- 
trix (IDM) of the generated ion. From the off-diagonal 
IDM elements, the relative phase and the degree of co- 
herence between the ionic states can be extracted, which 
are highly sensitive to the multichannel interactions oc- 
curring during the ionization process [44]. Attosecond 
transient absorption spectroscopy has also been used to 
study the dynamics of autoionizing states [45-48] and 
to study the motion of an electron wave packet during 
ionization [49]. 

The rapid technical advances in synthesizing light 
pulses made it possible to generate sub-cycle NIR pulses 
lasting no longer than a few femtoseconds, and to re- 
duce the jitter (time delay fluctuation) between the NIR 
pump pulse and the XUV probe pulse to tens of attosec- 
onds [35]. This time delay stability allows one to reli- 
ably probe the NIR-driven tunnel ionization dynamics 
within an optical cycle (approx. 2 fs) as a function of 
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the pump-probe delay. For non-overlapping pump and 
probe pulses, the transient absorption spectroscopy can 
be used to determine the instantaneous IDM at the time 
of the probe pulse [40] . 

The aim of this study is to investigate the ion popu- 
lation dynamics in krypton within the pump pulse. We 
show that the instantaneous ionic state population can be 
well captured by the transient absorption spectrum even 
for overlapping pump and probe pulses. Furthermore, 
we observe strong modifications of the absorption lines 
in the transient absorption spectroscopy when pump and 
probe pulses have the maximum overlap. 

We show that these deformations can be understood by 
relative phase shifts in the ionic dipole. Wc identify that 
the highly non-perturbative dressing of the A-electron 
states (particularly with the neutral ground state) is re- 
sponsible for the phase shift. Also the dressing of the 
ionic (A — 1-electron) states, which leads to energy shifts 
in the ionic states, contributes to the phase shift. The lat- 
ter one is, however, much weaker than the first dressing 
mechanism. Note that these two dressing mechanisms are 
quite different in nature. The first mechanism dresses N- 
electron states and the second one dresses N — 1-electron 
states. 

To capture these dressing mechanisms during ioniza- 
tion, a description of the entire A-electron system is re- 
quired. We describe the dynamics of the full iV-body 
wave function with a time-dependent configuration inter- 
action singles (TDCIS) approach [50]. The description of 
the pump and the probe steps of Ref. [35, 39] requires at 
least two active electrons, since the pump pulse ionizes 
an outer-valence electron and the probe pulse resonantly 
excites an inner-shell electron into the generated hole. 
Therefore, it is crucial to use a multi-channel model, 
which goes beyond the single-active-electron (SAE) ap- 
proximation, to describe the pump-probe process. 

The paper is structured as follows. In Sec. II, we give 
an overview of our TDCIS method, and describe the the- 
ory of attosecond transient absorption for overlapping 
pump and probe pulses. The results are discussed in 
detail in Sec. III. In Sec. IV we draw our conclusions. 
Atomic units [51] are employed throughout, unless oth- 
erwise noted. 

II. THEORETICAL METHODS 
A. Equations of Motion 

The time-dependent Schrodingcr equation of an A- 
electron system exposed to linearly polarized electric 
fields is given by 

i^m)) = !*(*)>, 

H(t) = Hu + Hi- E(t)z, (lb) 

where \^(t)) is the full A-electron wave function and 
H(t) is the exact A-body Hamiltonian, which can be 



partitioned into three main parts: (1) Ho = F — irjW is 
the sum of the time-independent Fock operator F and a 
complex absorbing potential (CAP), which reads W(r) = 
[r — r CAp] 2 ®(r — ?"cap)j where r is the radius, and 0(r) 
is the Heaviside step function; (2) the electron-electron 
interactions that cannot be captured by the mean-field 
potential in H a are captured by Hi (= Vc — Vhf — -Ehfj 
for a detailed description of these quantities see Refs. [50, 
52]); (3) the term E(t) z is the laser-matter interaction in 
the electric dipole approximation using the length form. 
The CAP within Hq prevents artificial reflections of the 
ionized photoelectron from the radial grid boundary and 
is located far away from the atom such that all processes 
close to the atom are unaffected by the CAP. This is 
controlled by the parameter tcap- 

By allowing only one electron to get excited or ion- 
ized out of the ground state configuration, we strongly 
reduce the complexity of solving Eq. (la). A suitable 
way to achieve this goal is by exploiting the configura- 
tion interaction (CI) language and describing the A-body 
wave function in terms of the Hartree-Fock ground state 
|$o) and singly-excited configurations |$f ). This approx- 
imation is known as Cl-Singles (CIS). The corresponding 
TDCIS A-electron wave function reads 

]*(*)) = ao(t)|*o> |*?>, (2a) 

i . a 

m) = ci^o), (2b) 

where i,j and a, b refer to occupied orbitals and unoccu- 
pied (virtual) orbitals, respectively, in the Hartree-Fock 
ground state |$o) • Indicies p, q are used when no distinc- 
tion is made between occupied and virtual orbitals. The 
operators c p and c p create and annihilate, respectively, 
an electron in the spin orbital \<p p ). The equations of 
motion (EOM) for the expansion coefficients a (t) and 
af(t) read: 

ia (t) = -£(t)£;(#o|W) <*?(*) (3a) 

i . a 

id?(t) = ( £a - £i )«?(i) + ^(^l^i|^)»iW ( 3b ) 
-£(t)(($?| z |$ ) a (t) + £ m z a){t)) , 

where e p are the orbital energies of the orbitals \<p p ), 
which are eigenstates of the modified, time-independent 
Fock operator, i.e., Hq \ip p ) = e p \<Pp). The operator Hi 
is the residual electron-electron interaction, which goes 
beyond the mean-field potential. The parentheses |-) 
and (-| stand for the vector and dual vector with re- 
spect to the symmetric inner product required because 
of the non-hermiticity of Hq. The dipole interaction be- 
tween singly-excited configurations reduces to transitions 
between states of the excited electron and transitions be- 
tween ionic states: 

(®t\*A®])=(v«\~ z \Vb)Ki-{Vj\Z\Vi)Kb- (4) 
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A detailed description of our implementation of the TD- 
CIS method can be found in Ref. [50, 53]. 

From the full iV-body wave function one can construct 
the ion density matrix (IDM) p lBM (t) by tracing over the 
excited electron. The matrix elements are given by 

pIT® =E («?(*) [«?(*)] *o b>a + 2r, e^^X (5) 

a,b 

x f dt' w bta af(t') [a){t')]* e- i{Ei ~ ei)t ' 

J — oo 

where Wb, a are the matrix elements of the CAP in the 
virtual orbital basis, and Ob, a are the overlap matrix el- 
ements between virtual orbitals. The second term in 
Eq. (5) corrects the loss of norm in the IDM due to the 
absorption of the excited electron by the CAP. The CAP 
is placed far away from the atom such that an electron 
so far out does not affect the ion, specifically the ionic 
states. Therefore, the absorption of an electron by the 
CAP results only in an artificial loss of norm that is com- 
pensated by the second term in Eq. (5). To keep the no- 
tation compact, we use the notation p] BM (t) := pYi M (t) 
for ionic state populations. 



the electron, and \n,l,m) are the 1-particle orbitals ob- 
tained from a non-rclativistic Hartree-Fock calculation. 
The orbital energies £j are taken from experimental ion- 
ization potentials. For the virtual orbitals, we can ne- 
glect spin-orbital splitting and use the quantum numbers 
na,la,Va, m>a to classify the orbitals, where a a is the spin 
component in the laser polarization direction and is 
the projection of the orbital angular momentum onto the 
laser polarization direction. The virtual orbitals read 



(7) 



The usage of linearly polarized light leads to the con- 
dition, mf = + a a , for each singly-excited configura- 
tion which is conserved by H(t). After the intro- 
duction of spin-orbit splitting for the occupied orbitals, 
we cannot make use of the a and m L symmetries inde- 
pendently to reduce the number of singly-excited con- 
figurations, |5>f), as done in Rcf. [53]. However, not 
all symmetries are lost and we find that Eq. (3b) is (up 
to a global phase) invariant under the parity transfor- 
mation (m(,m%,(T a ) —> (—mf,—m^,—a a ). The new 
parity-adapted, singly-excited configurations read 



B. Spin-Orbit Splitting 
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In order to include the effect of spin-orbit splitting 
in the occupied orbitals, wc follow the logic of Ref. 
[54], where we account for spin-orbit splitting with 
degenerate-state perturbation theory within the (n, I) or- 
bital manifold. The occupied orbital i is, then, charac- 
terized by the quantum numbers n,i,li,ji,mf , where rij 
is the principal quantum number, l{ is the orbital angu- 
lar momentum, ji is the total angular momentum, and 
mf is the projection of the total angular momentum onto 
the polarization direction of the external laser field. The 
occupied spin-orbitals read 



where the configurations |5>") ± stand for singly-excited 
configurations with mf ^ 0, respectively. States with 
7r = are gerade parity state and 7r = 1 states 
are ungerade parity state. All ungerade configurations 
|$") 7r=1 will never get populated and we can exclude 
them in our further investigations. Note that the fac- 
tor ( — 1 )'i+ s « — ii comes from the symmetry Cf 3 'J™ 3 , „, = 

( — \\li+h—h f<h,—m 3 

V / li ,— fni ;l2 , — rn,2 ' 



li f mf + i ;si 



l \l^i> hi mf 2) 

1 \n,i,li,mf + 5) 



(6) 



where the Clebsch-Gordan coefficient is given by 



C 3 
1 1 



(h, mi; I2, ?7i,2^3, ^3), Si = 5 is the spin of 



Matrix Elements 



The matrix elements, which are needed for Eqs. (3), 
must be evaluated in the parity-adapted \&f) n configu- 
ration basis. The dipole matrix elements with respect to 
the symmetric inner product read 



^ ^7Tl ,7T2 •) 

(*o| = ^,0*m/,mJ+«r.V^«^j ) C^ > 5_ (r(t . a4iffB , (9b) 

where we used the notation i = mf) with I = (n,, U) and a — (A, m„, cr a ) with A = (n a , l a ). The dipole matrix 

elements in the original, non-parity-adapted basis are given by z™ p := (n p , l p , m\ z \n p , l q , m). The matrix elements 
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of Hi read 



= 2 5 nifi 5 T2t0 d i ' m i c 3 J' m i v^ nm (10) 



\AJBI) 



where Mi = (m^m^m^mf;), Mf = (m^mf - a, m^,m{ - cr), and = (m^-mf - a, -m£, m/ - a). The 
Coulomb matrix elements in the non-parity-adapted basis read v^Q RS y=(n p ,l pi m p ;n q ,l qi rn q \ri 2 \n r ,l r ,m r ;n Sl l s ,m s ) 
with M = (m pi m q ,m r ,m s ). 



D. Transient Absorption for Overlapping Pulses 

The transient absorption signal is a direct measure of 
the cross section of the system. In Ref. [40] the tran- 
sient absorption signal was derived for the case of non- 
overlapping pump and probe pulses. The probe pulse 
was treated in first-order perturbation theory such that 
it was possible to give an analytic expression for the tran- 
sient absorption signal as a function of the instantaneous 
IDM p lBM (t). The pump pulse, usually a strong- field 
NIR pulse, which ionizes the atom by tunnel ionization, 
was treated non-perturbatively. 

For overlapping pump and probe pulses, the influence 
of the probe pulse does not decouple from the impact of 
the pump pulse. Therefore, it is not clear to which extent 
p IDM (i) can be extracted from the transition absorption 
spectrum like for non-overlapping pulses. In order to 
fully capture the overall effect of pump and probe pulses, 
both pulses are treated non-perturbatively meaning the 
TDCIS equations of motion [see Eq. (3)] are solved for an 
electric field E(t) = E pump (t) + E pm b e (t). Note that the 
probe step could also be treated pcrturbativcly by intro- 
ducing a two-time IDM, which depends on two different 
time arguments. In our non-perturbative approach only 
the one-time IDM p lBM (t) needs to be constructed for 
each pump-probe configuration. From p IDM (i) the ionic 
dipole moment 



IDM, 



and the atomic cross section 



a a {uj) = 4-7T — Im 
c 



(*)io,M 

-Eprobe(w) 



(11) 



(12) 



can be calculated. By performing the trace over /5 IDM (i) 
and not over the full iV-body density matrix p(t), we con- 
sider only dipole transitions between ionic states. Transi- 
tions between virtual orbitals can be neglected, since the 
XUV probe pulse interacts only weakly with the excited 
electron. Transitions between occupied and virtual or- 
bitals describe stimulated emission and photoionization 
processes. Both mechanism do not lead to sharp fea- 
tures in a a {u) around the bound-bound transition en- 



ergies. Therefore, we ignore these contributions, which 
lead to background signals we are not interested in. 

The detector, where the transient absorption spectrum 
is measured, does not record the atomic response but 
rather a damped spectrum of the form 



\E pro bo(L 7 w)| = |-Ep ro be(0, W)| 



(13) 



where -E plo bc(0, w) is the incoming probe electric field, L 
is the length of the medium, E pro b c (L, uj) is the probe 
electric field at the end of the medium, and tiat is the 
atomic number density. In Eq. (13) Beer's law is used, 
which assumes a homogeneous medium and that the ra- 
tio ( z )ion M/^probcM is independent of E pmhc (cu) := 
£'probe(0, w). In Sec. IIIB the validity of Beer's law is 
discussed. 

Due to the finite energy resolution of the detector, 
the transient absorption signal in Eq. (13) has to be 
convolved with a Gaussian mask function, where the 
full-width-at-half-maximum (FWHM) width is given by 
the energy resolution of the detector. The cross section 
(7 m (w) measured at the detector can be related to the 
atomic cross section and is given by 

<r m {u) = l —\n ( e -»AT W«) * G SE (ojj) , (14) 

71 AT L \ J 

where Gse^) is the area- normalized Gaussian with the 
FWHM width of SE, and the symbol * stands for the 
frequency convolution. The dependence on the pump- 
probe configuration, specifically the pump-probe delay 
r, enters parametrically in Eqs. (11)-(14) such that the 
atomic and the measured cross sections read a a (u>; r) and 
cr m (w;r), respectively. 



E. Oscillating Dipole Model 

In the following, we develop general expression for the 
transient absorption spectrum, which is based on a simple 
model. Later in Sec. Ill we use this generalized expres- 
sion to discuss the features of the transient absorption 
spectrum for overlapping pulses obtained from our TD- 
CIS calculations described in Sec. IIA-IIC. 



5 



First, we reduce the description of the ion to a two- 
level system. The ground state \g) can only be ac- 
cessed by the pump pulse via tunnel ionization and the 
excited state |e) can only be accessed by the probe 
pulse via resonant excitation out of \g). The probe 
pulse, which may be approximated by a delta pulse, 
i.e., E VTO \, e [t\T) = E()5(t — r), creates a coherent su- 
perposition > t)) = a Q \g) + ai e -*(<"o-ir/2) (t-r) 
where uq is the positive energy difference between the 
two states, 1/r is the lifetime of the excited state, and 
ai = —iE a (e| z \g) results from the excitation by the 
probe pulse. This superposition leads to an oscillating 
dipolc 

<^>io„ (* > t) = m)\ z Mt)) = -2E a | (e\ z \g)\ 2 \a \ 2 
x sm[uj Q (t -T)}e-^^- T \ (15) 

Inserting Eq. (15) in Eq. (12) and using E pTO h e (<jJ', t) = 
EQe~ tulT the final expression for the cross section reads 



o-(w;r) 



All U! 



A) 



r/2 



(w - w ) 2 + V 



(16) 



where zo = l a o| 2 |(e| z determines the transition 
strength. We see that for a simple two-level system the 
cross section is purely Lorcntzian and directly propor- 
tional to the ground state population |ao| 2 at the time of 
the probe step. 

Adiabatic energy shifts in the ionic states during the 
intense NIR pulse result in a phase shift in the oscillating 
ionic dipole, i.e., (z) ion oc sin[cjo(^ — t) + (/>(t)]. Here, we 
assume the ionic state and the dipole oscillation live for a 
long time after the NIR pulse is over such that the entire 
dipole dynamics can be approximated by a phase-shifted 
oscillation. The phase shift 4>{t) has a dramatic influence 
on the shape of the transition line, which reads 



<j(uj;t) 



4ttlo t cos[0(t)] + (uj - ujq) sin[</>(r)] 
z - — 



(w - ujo) 2 + V 



(17) 



Note that the phase shift 4>(t) affects only the shape of 
the transition but not the strength zq. 

In Fig. 1 cr(oj;r) is shown for specific values of <fi. The 
transition line is purely Lorentzian for <f> = 0. In the 
case 4> = 7r/2, the cross section shows a dispersive behav- 
ior and has equally negative and positive regions that 
lie symmetrically around the field-free transition energy 
wo(0). For all other phases, the cross section is a sum of 
these two scenarios and becomes asymmetric around wq. 
A phase shift by it changes the sign of the cross section. 
For — 7r/2 < <fi < 7r/2, the system shows an absorbing be- 
havior whereas for7r/2<0<3/27r the system is rather 
emitting. A similar scenario has been discussed in atomic 
helium, where neutral excited states are dressed by an IR 
pulse leading to emitting and absorbing patterns depend- 
ing on the pump-probe delay [55]. 
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FIG. 1: (color online) The cross section a(tj) of Eq. (17) is 
shown for several values of 4>(j). The transition energy is 



ljo = 3 a.u., r = 3.2 x 10" 
is given by zq = 0.01. 



a.u., and the transition strength 



Similarly to the dressing of the ionic states, the in- 
fluence of the excited electron on the ion via the resid- 
ual Coulomb interaction and via the pump field can 
lead to additional phase shifts in the oscillating dipole 
(see Sec. II F). Furthermore, corrections to the transition 
strength zq can occur, which may cause zq to be no longer 
directly proportional to the instantaneous hole popula- 
tion. In order to capture these effects, we parametrize in 
addition to the phase 0(r) also the transition strengths 
zq — > Zq{t). Our generalized version of Eq. (17) for a 
multi-level ion reads 



Cdipole(w; t) 



47TW ^ 

> z t (t) 

c ^— ' 

T 



(18) 



cos[(/> t (t)] + (u - uj t ) sm[(j)T{r)] 

where we sum over all possible ionic transitions T. Note 
that Eq. (18) is designed to capture the influence of all 
these effects that go beyond our simple two-level model 
[see Eq. (16)]. However, Eq. (18) cannot explain why 
these changes occur and where they come from. 



F. Mechanisms leading to the Phase Shift 

As discussed in Sec. HE, the dressing of the ion can 
induce a phase shift in the ionic dipole. In the following, 
we discuss in the language of TDCIS how the dressing 
by the field and the coupling of the excited electron to 
the ionic subsystem can influence the phases 4>t{t) of 
Eq. (18). First, we analyze the scenario where the time 
evolutions of the excited electron (a index) and the ionic 
states (i index) are decoupled. This is the case when 
the terms ($£]#! and (*il^l*o) 

are switched off 
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in Eq. (3b). The resulting EOM can be written as 
ia« (t) = ^ Hf:* b) (t) a\ (t) + Aft}) (*) a l (*) > 



where H? lc ?, 

(a,b) 



(t) := e b S a . b +E(t) (cp a \z \ip b ) affects only the 

excited electron, and flS> (i) := — 5^j—E{t) (ipj \z\ipi) 

affects only the ionic states. H Qlcc and H ion can be 
viewed as Hamiltonians of the two subsystems. Note 
that the Hamiltonians of both subsystems commute, i.e., 

H elec (t), H ion (t) = 0. To confirm that Eq. (19) leads to 
decoupled EOMs for the excited electron and the ion, we 
make the product ansatz af(t) = Xa(t)ni(t), where the 
EOMs of the separated electronic and ionic wave func- 
tions are given by 

6 

iki(t) = 2 #{$)(*)«*(*)■ ( 20b ) 

3 

We find that the product ansatz with the decoupled 
EOMs of Eq. (20) solves Eq. (19). This shows that 
Eq. (19) is the overall EOM of the full TV-electron sys- 
tem, which consists of two totally separated subsystems. 
Note that the term [cf. Eq. (4)] does affect 

the excited electron and the ion but it does not lead to 
interactions between the two. Enforcing a normalized 
electron wave function, we find from Eq. (20b) that the 
IDM is given by 

(21) 

where T = i — >• j denotes the ionic transition. If E(t) = 
for all times t, we find 4>^ n (t,T) = for all t and r. 
Hence, the additional phase 4>^ n (t, r) enters only due to 
the term —E(t) (<fj \ z \ipi) in H lon . This is exactly the 
field-driven dressing of the ionic system. After the pulse 
is over [E(t) = 0], the phase 0^ n (i, r) becomes indepen- 
dent of t and depends only on the probe time r, i.e., 
0jp n (t,r) -> 0j? n (r). In Sec. IIID1, we show how the 
field-driven dressing of the ionic system can be analyti- 
cally analyzed with the help of the polarizability of the 
ion. It is interesting to note that the influence of the 
electric field on the excited electron, i.e., E{t) (<p a \ z \<fib), 
does not influence the IDM and subsequently the ionic 
dipolc oscillation. 

Additional phase shifts similar to 0^ n (r) can also oc- 
cur due to the coupling between the ion and the excited 
electron. There exist two kinds of mechanism that can 
couple these two subsystems: (1) the residual Coulomb 
interaction [the ($°|7?i in Eq. (3b) are the corre- 
sponding matrix elements], and (2) the field-driven mix- 
ing of the excited ^-electron states with the neutral 
ground state [the ($"|£|$q) m Eq. (3b) are the corre- 
sponding matrix elements]. Both terms were ignored in 



Eq. (19), which led to two decoupled subsystems. To 
distinguish the phase shifts induced by the two differ- 
ent mechanisms, we introduce 4>t sldual (t) and </)|f oun (r). 
The phase shift due to the residual Coulomb interac- 
tion is denoted by </>5? sldual (r), and i^ff oun (r) denotes the 
phase shift due to the field-driven mixing to the neutral 
ground state. 

Adding up all three phase shifts we find that the total 
phase shift for the transition T is given by 



Mr) = ^"(r) + ^ csidual (r) + ^ round (r). 



(22) 



In Sec. HID we discuss which of the three phase shifts 
gives the dominant contribution. 



III. RESULTS 
A. System, Pulse, and Numerical Parameters 

All presented results were calculated with the XCID 
package [56]. All calculations were performed with 600 
radial grid points, a maximum radius r max = 150, and 
a non-linear grid mapping parameter £ = 1.0. The 
CAP starts at a radius tcap = 130 and has a strength 
r\ = 0.003. The maximum angular momentum employed 
was i max = 30. Furthermore, if the orbital angular 
momentum of any 1-particle orbital \(p p ) appearing in 
the matrix elements of H± exceeded 4, then only the 
monopolc term of H± was considered. With these param- 
eters, our results are converged. A detailed explanation 
about the parameters of the grid, of the CAP, and of the 
residual Coulomb interaction can be found in Ref. [50]. 

The NIR pulse profile measured in Ref. [35] is used as 
the pump pulse in all the following results and is shown in 
Fig. 2. The pump pulse is approx. 2 fs long with one main 
and two side peaks. The maximum electric field strength 
of the pump pulse is ps 0.08 a.u. (corresponding to an 
instantaneous intensity of 4.8 x 10 14 W/cm 2 ). For the 
probe pulse (also shown in Fig. 2) a Gaussian pulse profile 
is used with a central frequency of 3 a.u. (w 81 eV), a 
FWHM- width of the intensity profile of 10 a.u.(« 240 as), 
a carrier envelope phase of zero, and a maximum field 
strength of 10~ 2 a.u. (ra 3.4 TW/cm 2 ). 

As in Refs. [35, 39], we choose atomic krypton as our 
system of interest. The spin-orbit coupling within the 
occupied orbitals is accounted for in first-order perturba- 
tion theory as described in Sec. II B. The energies of the 
singly- ionized ionic states [NLj 1 ]^ 1 with respect to the 
Hartree-Fock ground state energy are given by the neg- 
ative orbital energies (i.e., Koopmans' theorem). Here, 
N is the principal quantum number, L the orbital angu- 
lar momentum, J the total angular momentum, and M 
the projection of J. Since we are using linearly polarized 
light, the values for M and — M states are always the 
same. Therefore, we refer always to a sum of both \M\ 



contributions, i.e., \NLf\~ 1 = [NL^]' 1 + [NT 
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FIG. 2: (color online) NIR pump (red solid) and XUV probe 
pulses (green dashed) used for transient absorption spec- 
troscopy are shown separately. The peaks of both pulses are 
centered at t = 0. 



The orbital energies of the Ap orbitals are replaced 
by experimental values taken from Ref. [57] in order 
to match the experimental ionization potentials. The 
orbital energies of the 4p-shell are given by £4 Pl/2 = 
-0.5389 (= -14.67 eV) and e 4p3/2 = -0.5148 (= 
— 14.00 eV). The orbital energies of the 3d-shell are 
taken from Ref. [58]. In addition, we account for the 
finite lifetime of all 3d -1 configurations of 7.5 fs (T = 
88 meV) [39]. Hence, the 3d orbital energies become 
complex and read £3d 3/2 = —3.5 — iT/2 (= —95.24 eV — 
i0.044 eV) and e 3d&/2 = -3.4525 - iT/2 (= -93.95 eV - 
i0.044 eV). 



B. Transient Absorption Spectrum 

In a recent experiment [35] attosecond transient ab- 
sorption spectroscopy has been used to investigate the 
hole production dynamics in atomic krypton during a 
sub-cycle NIR pump pulse. The theory described in 
Ref. [40] was used to analyze the transient absorption 
spectrum and to connect the spectrum to the instanta- 
neous IDM pYj M (t). Propagation effects and the finite 
energy resolution of the detector were accounted for as 
described in Eq. (14). The macroscopic propagation is 
captured by Beer's law. Previous studies [40] have shown 
that Beer's law is valid for the pump-probe scenario in- 
vestigated here. Similar conclusions have been found in 
Ref. [41] when the probe pulse is much shorter than the 
pump pulse (as in our case). For probe pulses longer than 
the pump pulse, the macroscopic propagation can quite 
strongly deviate from Beer's law. 

The theory developed in Ref. [40] is strictly speaking 
only valid when pump and probe pulses do not overlap 
and the tunnel-ionized electron is far away from the par- 



ent ion. In this case, the ionic subsystem, i.e., p lBM (t), 
reduces to a simple multi-level system without any kind 
of interactions with the environment and between lev- 
els. Hence, the entire dynamics of the ion is analytically 
known and reads p^ u (t) = /4° M (io)e i(£i - E ' )( *-* o) . For 
overlapping pulses, the dynamics of p lBM (t) become more 
complex. Therefore, it is not clear to which extent the 
transient absorption spectrum for a given pump-probe 
delay can be related to pYj M {t) when pump and probe 
pulses do overlap. 

In Fig. 3 the calculated transient absorption spectrum 
cr m (w;r) is shown as a function of photon energy u and 
pump-probe delay r. The three main transition lines, 
i.e., 4p~L -> Zd~j 2 (79.95 eV), 4p~l -> 3d~l (80.57 eV), 

and 4pg/ 2 — > 3dg, 2 (81.24 eV), are clearly visible. To 
shorten the notation we refer to these three transition 
lines as Ti,Ts, and T3 , respectively. 

The cross section shown in Fig. 3 is in agreement with 
experimental observations [35]. The transition strengths 
increase mainly around r « 0, when the krypton atom is 
being probed within the main peak of the pump pulse. 
It is also during this main peak of the pump pulse where 
the atom gets mainly ionized. Simultaneously during the 
hole creation, the transition lines in the transient absorp- 
tion spectrum change dramatically their shapes, resulting 
in negative cross sections for energies just below the field- 
free transition energies. Negative cross sections can be 
only seen here for the T\ transition, because a detector 
resolution of ~ 300 meV let these features disappear for 
the transitions T 2 and T 3 (cf. Fig. 4). 

Similar to r m 0, line deformations and negative cross 
sections do also occur at to ss 80 eV for r » ±1 fs, where 
the probe pulse coincide with the side peaks of the pump 
pulse (cf. Fig. 2). The side peaks of the pump pulse are 
strong enough to lead to tunnel ionization as well (cf. 
Fig. 5). Particularly the ionization and the deformation 
caused by the first side peak of the pump pulse (r 
— 1 fs) can be nicely seen in Fig. 3 for the transition line 
T\ . As it will become clear in the discussion in Sec. Ill D, 
the mechanism behind the deformations in all three cases 
(r w 0, ±1 fs) is the same. 

For t > 2 fs, the dynamics of the hole populations 
are barely affected by the pump pulse and behave as if 
they were field-free. At these time delays, the intensity 
of the pump pulse is also strongly reduced (by more than 
a factor 10) compared to the peak intensity at t = 0, 
thus supporting the observation of field-free behavior for 
larger pump-probe delays. Field-free behavior means for 
the main transition line T\ that it becomes stationary 
and does not change anymore in shape and strength. The 
other two transition lines show interference effects from 
the coherent superposition of Ap~j 2 and Ap~J 2 [39, 40]. 

Large negative time delays (r < — 2 fs) are not of in- 
terest and, therefore, they are not shown in Fig. 3 due to 
two reasons. First, the ionic cross section [see Eq. (12) 
and Eq. (14)] is zero for large negative time delays, be- 
cause only neutral krypton atoms exist prior to the NIR 
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FIG. 3: (color online) Attosecond XUV transient absorption spectrum a m (ijj;r) [see Eq. (14)] of krypton as a function of 
energy ui and pump probe delay r. 



r = 70 (~ 1.7 fs). The transition lines in a m (uj;T) 
are broadened with respect to the atomic cross section 
<7 a (cj;T) due to the propagation effect and the finite de- 
tector resolution (sa 300 meV), which is wider than the 
natural transition widths (I\ = 88 meV). The cross 
section <7dipoic(^; t) is obtained by fitting Eq. (18) to 
cr a (uj; t) obtained from the TDCIS calculations. The en- 
ergies ujTi and the natural widths Tt { of all transition 
lines Tj are kept fixed (see Sec. Ill A) and only the magni- 
tudes ZTi{r) and the phases ^t^t) are fitted to er a (w; r). 
The features of a a (ui; r) are well captured by <7dipoic(w; r) 
for all t. At r = 70, the two strongest transition lines are 
Lorentzian-shaped as expected for non-overlapping pump 
and probe pulses. The second (Xb) and especially the 
third (T3) transition lines do not have a Lorentzian shape 
due to the coherent superposition of 4p~^ 2 and 4p7/ 2 . 

The success of <Jdipoio(w; r) in capturing all features of 
a a (oj; t) shows that the influence of all terms in Eq. (3), 
which go beyond a simple two-level model (see Sec. II F), 
can be understood by phase shifts </>T; (t) and changes in 
the oscillating dipole strengths (t) . 

However, the oscillating dipole model cannot explain 
what is the physical origin of <f>Ti (j) or whether or not 
z Ti(j) can be related to p IDM (r). The answers to these 
questions are discussed in the following. 



pulse. Second, no electronic dynamics can be probed be- 
fore the NIR pulse, since the neutral krypton atoms are 
in the electronic ground state. 
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photon energy w [eV] 
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FIG. 4: (color online) Transient absorption spectra for the 
pump-probe delay r = and r = 70 (f» 1.7 fs) are shown 
in panel (a) and (b), respectively. The atomic cross section 
<j a (red solid), the measured cross section a m (green dashed), 
and (Tdipoio (blue dotted) are shown for each r. 



In Fig. 4 the transient absorption spectra cr a (oj;r), 
cr m (w;r), and 

^"dipole (w,t) are shown for r = and 
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C. Population dynamics 

First, wc turn our focus to the population dynam- 
ics of the ionic states. Particularly, we investigate the 
hole creation dynamics in the 4p 3 / 2 orbitals during the 
pump pulse. In order to do so, we need to focus only 
on the main transition T\. For non-overlapping pulses, 
the transition strength is proportional to the instan- 
taneous population p™ M (r) := pf™ 2 { T ) + fp^(r), 

where /0* dm (t) are the hole populations of the ionic states 
i = [4p^2] 1 and i = [Apy^] _1 (for details on the no- 
tation see Sec. Ill A), respectively. In Fig. 5, we com- 
pare pJ,g M (r) (red solid) with the transition strength 
^drpoic ,_ ZTi ( T } (blue dotted), with the reconstructed 
population p™n (green dashed) obtained by applying the 
same fitting procedure as in Ref. [35] to rr m (u;;T), and 
with the experimental data of Ref. [35] . 
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FIG. 5: (color online) The instantaneous hole population 
Pcff M ( r ) ( re d solid) is shown together with the reconstructed 
populations Pcff P ° 1C ( r ) (blue dotted) and p^g (green dashed). 
The experimental population (grey dashed) is taken from 
Ref. [35]. The reconstructed and experimental populations 
are scaled such that all have the same value at t = 2.4 fs. The 
NIR pulse intensity is highlighted (gray area) in the back- 
ground. 



In all three theoreticcal curves, oscillatory behavior can 
be seen during build up of the hole populations. The 
features of p|,g M (r) are well captured by p^ p ° 1c (t) and 
by Pch(t). It shows that p 1 J^ m (t) can be quite well re- 
constructed from the transient absorption spectrum even 
though pump and probe pulses do strongly overlap. Note, 
however, that at the major peak (r £=* 0) a delay of about 
w 200 as occurs in p^(r). The deviations from pJ.§ M (r) 
are a measure of the influence of the excited electron on 
the ionic states, since the adiabatic dressing of the ionic 
states does not affect the transition strength. As de- 
scribed in Sec. II E, if the evolution of the excited electron 
decouples directly after the probe step from the evolution 



of the ion, p^ pole (r) coincides with pJ.§ M (r). 

In the experiment, this oscillatory behavior was not 
seen and the transition strengths increased monotonically 
with r. The reason for this discrepancy might lie in the 
strong restriction of the CIS space for the ionic degree 
of freedom. Within the CIS space generated from the 
neutral ground state of the atom, the ion is described 
exclusively by one-hole configurations. 



D. Line Deformations and Phase Shifts 

As we have seen in Sec. IIIB the oscillating dipole 
model is able to describe all features of the transient ab- 
sorption spectra by only adjusting the strengths Zt(t) 
and the phases <Pt(t) of the three transitions. The 
strengths Zt(t) are in close connection to the instanta- 
neous populations of the ionic states even for overlapping 
pulses [cf. Sec. IIIC]. The correct values of <Pt(t) arc 
important to capture the shapes of the transition lines, 
which can change significantly during the pump pulse [cf. 
Fig. 4]. 




-10 12 
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FIG. 6: (color online) The phase shift 4>t^ (t) of the strongest 
transition line and the phase shift </>t 3 (t) of the weakest tran- 
sition line are shown. 

In Fig. 6, </>Ti( r ) an d 4>t 3 ( t ) arc shown. The phase 
0t(t) is chosen such that the strengths zt are always 
positive. The region r < —1.5 fs is not shown, since 
the strengths of the transitions are of the same order 
as the numerical background noise, which leads to large 
uncertainties in </>j\ The last peak of the pump field at 
t = 1 fs has only a small effect on 4>t whereas the first 
peak at r = — I fs has a strong influence on the phases. 
This is an interesting observation, since both peaks have 
the same intensity. 

As mentioned earlier, the coherent superposition of 
the ionic states Ap^, 2 and Ap~J 2 results in a phase shift 
( £ 4p3/2 — £ 4pi/i) T i which is particularly dominant in tran- 
sition T3 [40] and has been observed in previous cxper- 
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iments [39]. This coherence-induced phase shift is also 
visible in 4>t 3 ( t ) for t > 1 fs (see Fig. 6). The slope of 
4>t 3 (t) is approx. 0.023 a.u. which is in good agreement 
with £4p 3/2 — £4p 1/2 = 0.024. The main transition line T\ 
is not affected by this coherent superposition and, there- 
fore, goes over into a Lorentzian shape (i.e., 0Ti = 0) for 
large r (cf. Fig. 6). 

We have seen in Fig. 6 that the phase shift 4>t can be 
quite large. However, what has not been answered yet is 
the origin of 4>t- As discussed in Sec. II F, there are three 
main contributions to 4>t{t) [cf. Eq. (22)]: (1) </>y n (r) in- 
duced by the NIR-driven dressing of the ionic system [the 
(ipj\z\(pi) are the corresponding matrix elements]; (2) 
(/>55 sldual (r) induced by the residual Coulomb interaction; 
and (3) 0|, lound (r) induced by the NIR-driven mixing of 
the excited TV-electron system with the neutral ground 
state [the ($°| z |$o) are the corresponding matrix ele- 
ments]. Particularly, in order to account for 0^ sldual (r) 
and </>| n round (r) , it is important to have a multi-electron 
picture, which can describe the degrees of freedom of the 
ionized electron and of the ion. 
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FIG. 7: (color online) The apparent energy shifts of the 
strongest transition line Tj are shown for the calculated cross 
sections with (red line with crosses) and without (green with 
asterisks) the 4s orbital active, and for the experimentally ob- 
tained cross sections (blue dotted). The experimental values 
are taken from Ref. [35]. 



The strong line deformations, if only the positive part 
of the cross section is considered, appear as if the tran- 
sition energies have moved. This energy shift we call ap- 
parent energy shift. In Fig. 7, we make a direct compar- 
ison between calculated (red and green lines) and exper- 
imentally obtained (blue dashed line) [35] energy shifts 
of the strongest transition line, T\. In the calculations 
yielding the red line (with crosses), no approximation is 
made in the TDCIS calculations. For the green line (with 
asterisks), all couplings to the 4s orbital are turned off. 
The procedure used to extract the apparent energy shift 
is described in Ref. [35]. 



The magnitude of the apparent energy shifts are cor- 
rectly reproduced by the calculations. Ignoring the 4s 
orbitals makes no significant difference. Particularly in 
the theoretical results, the second pump peak is clearly 
visible. For r > 1.4 fs, the energy shifts go as expected 
to zero indicating that influence of the pulse decreases 
and the krypton ion can be treated as a field-free ion. 
Negative pump-probe delays are not shown because the 
transition strengths are too weak to obtain reliable re- 
sults. 

In the following, we discuss each mechanism which can 
lead to the strong line deformations seen in Fig. 3. We 
discuss each mechanism in terms of 4>t(t), since the oscil- 
lating dipole model showed that these deformations can 
be fully understood by a phase shift. 

1. Field Dressing of the Ionic System and the Polarizability 

of Kr + 

In the presence of a strong external field, the energies 
of the ionic state become modified, i.e., Si — > £i(E). If the 
photon energy is much smaller than the ionization poten- 
tial, the influence of the field on the ionic state can be 
described adiabatically. This is the case for Kr + , where 
the ionization potential is 24 eV [59], and the photon 
energy of the pump pulse is 1.4 cV and, therefore, far 
off-resonance with any ionic transition. Hence, the en- 
ergy corrections Aei(E) = Si(E) — £j(0) of the dressed 
ionic states i are well captured by the quadratic Stark 
effect 

Ae, t [E„(t)} = ^ £ p 2 ump (i), (23) 

where a(i) is the polarizability of the ionic state i, and 
the pump electric field is linearly polarized. Note that 
the energy of the ionic state i is — After the pulse 
is over, Aej(0) = and the ion is back in its field- free 
state with the ionic energies £;(0). If an oscillating dipole 
(coherent superposition of two states i and j) is present 
while the energy of the ionic states get shifted, the dipole 
oscillates after the pulse with the field-free transition en- 
ergy Wj_> 3 (0) = £j(0) — £j(0) but phase-shifted by [35] 

/oc 
dt' A^[E punip (t% (24a) 

Auii^iE) ^Ae^E) - Ae 3 {E) = ^'^^ E\ 

(24b) 

where Awj_>.j[.Ep Um p(f)] is the instantaneous shift in the 
transition energy between the states i and j. The phase 
shift <j) l °^^{t) is exactly the ionic phase shift </>y n (r) with 
T = i — > j, which originates from the field-dressing of 
teh ionic system. 

In order to estimate <j)]f, n {i) we use the quantum chem- 
istry code DALTON [60] to calculate exact polarizabili- 
ties of the ionic states. Furthermore, we can clarify 



11 



TABLE I: Static dipole polarizabilities a x<x = a y ,y and a z ,z 
of several states of Kr + are shown. Polarizabilities obtained 
by the CIS and CASSCF methods are compared. The CIS 
and CASSCF calculations are done with dalton. All values 
are given in atomic units and with a precision up to the second 
digit. 







CIS 




CASSCF 






OLx ,x 


Otz,z 


OLx ,x 


Otz,z 




-i 


1.57 


0.01 


10.62 


10.77 




-i 


0.53 


2.09 


10.72 


10.57 




-i 


1.05 


1.05 


10.67 


10.67 




-i 


-0.01 


0.00 


9.57 


9.71 




-i 


0.00 


0.00 


9.63 


9.60 




-i 


0.00 


-0.01 


9.65 


9.54 




-i 


0.00 


0.00 


9.58 


9.68 




-i 


0.00 


-0.01 


9.65 


9.55 



to which extent our TDCIS calculation can correctly 
describe the polarizability of the ion and subsequently 
<Pi^lj(t). The polarizabilities were calculated with a com- 
plete active space self-consistent field (CASSCF) wave 
function and with a CIS wave function. Here, CIS means 
that Kr + is described in the space of one-hole configura- 
tions, which is consistent with out TDCIS configuration 
space. In CASSCF calculations, spatial deformations of 
the ionic states are included that would require CISD 
and higher order configuration excitations with respect 
to the ground configuration of the neutral atom. These 
deformations are not included in CIS. 

In Table I, we summarize the results of the static dipole 
polarizabilities for several states of Kr + with holes in 
the 4p or in the 3d orbital manifolds. For more details 
about the polarizability calculations see Appendix A. 
The CASSCF calculations are in good agreement (±2% 
accuracy) with the static polarizabilities in Ref. [61]. 
Polarizabilities obtained with the CASSCF method are 
in very good agreement (< 2%) with experimental re- 
sults [62, 63] for neutral krypton atoms. For ionic kryp- 
ton the correlation effects are reduced in comparison to 
neutral krypton, thus making the CASSCF calculations 
even more accurate. Hence, we may assume that the 
CASSCF results for Kr + are practically exact. 

The polarizabilities obtained with CASSCF have val- 
ues around 10.7 for ^pj 1 ionic states and values around 
9.6 for 3d" 1 ionic states. In both cases, the anisotropy 
(ctz,2 — o-x,x) is small. This stands in contrast to the CIS 
results, where a high anisotropy in the polarizabilities of 
4p is found. Furthermore, the polarizabilities of the 
ionic states 3d -1 are basically zero with the CIS basis 
set. The only contribution to a(3d ) comes from the 
weak coupling between the 3d and 4p orbitals. Within 
CIS, the polarizabilities for 4p _1 are determined by the 
coupling to the 4s _1 state. Ionic states with two or more 
holes in the initially occupied orbitals do not exist in CIS 



TABLE II: Differences between the static dipole polarizabil- 
ities a z ,z of the states involved in the main transition line 
~~ * 3^5/2- The values are given in atomic units. 



"(Wa]" 1 ) - «([3C/2] _1 ) 


CIS 


CASSCF 


l 

m = ~ 
m= § 


2.10 


1.03 


0.01 


1.17 



and, therefore, cannot contribute to the polarizabilities. 
This restriction in the CIS space limits the ability of the 
ionic states to respond to the external field and leads to 
much smaller polarizabilities (as seen in Table I). 

The quantity that determines <^ n (t) is the difference 
in the polarizabilities, not the polarizabilities themselves 
[see Eq. (24)]. We are only interested in the a zz com- 
ponent, since we use light linearly polarized along the z 
axis. For the strongest transition line 4p~^ 2 — > 3d~J 2 , we 

need to look at the differences a([Apy 2 ]~ 1 )~a([3d 7 ^ 2 ]~ 1 ) 
for the ionic states with m = 1/2 and m = 3/2. 

In Table II theses differences are shown. The strong 
m-dependence of the CIS results is due to the high 
anisotropy of the polarizabilities. The accurate CASSCF 
results show almost no dependence on m. Since the 
[4p^,~ 3 ^ 2 ] -1 population is in our calculation much smaller 

than the [^Py 2 ] _1 population, we focus only on the 
rn = | results. Contrary to the results of the polarizabil- 
ities, where the CIS approach underestimates the values, 
the difference a^p^] -1 )— a ([3d5/2] _1 ) i s overestimated 
by CIS. The maximum energy shifts in the transition for 
the given pump parameters [cf. Fig. 2] are 183 meV 
(CIS) and 102 meV (CASSCF), respectively. Assuming 
the maximum energy shifts persist over 2 fs, the resulting 
phase shifts are 0.56 rad (CIS) and 0.31 rad (CASSCF), 
respectively. This approximation can be nicely verified 
by a 3-levcl Bloch model describing only the (N — 1)- 
electron ionic states 3d - 1 ,4s _1 , and 4p _1 . 

The phase shift is, however, much smaller than 
(pTx [see Fig. 8]. Hence, <fr™(t) cannot be the main con- 
tribution to 0Ti(f)- Furthermore, when the coupling to 
the 4s _1 ionic state is switched off and a(i) ~ within 
TDCIS the phase shift <\>r x (t) is almost unchanged. Note 
that in this case, the ionic states cannot be dressed and 

= o. 

2. Residual Coulomb interaction 

Here, we discuss the influence of the residual Coulomb 
interaction — specifically (/)^5 sldual (r) . First, we simplify 
the krypton atom to a two-level system with the states 
[4P3/2] 1 an d [3^5/2] 1 such that no dressing of the ionic 
states can occur (</>^ n = 0). The total phase shift reads 
now t (t) = 05? sidual (r) + 0| round (T). Second, we also 
ignore Hi such that the ionized electron can have no in- 
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pump-probe delay t [fs] 

FIG. 8: (color online) The phase shifts 4>T y (t) of the transition 
line Ti(4pr, 2 — s> 3d7, 2 ) are shown for TDCIS calculations with 

(red solid) and without (greeen dashed) the ionic state 4s . 
Additionally, the phase shift $f? n (r) (blue dotted) is shown for 
a Bloch model describing 3 ionic states (see text for details). 
Here, (^"(t) originates solely from the dressing of the ionic 
state [4pgy2] _1 through the coupling to 4s _1 . 
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FIG. 9: (color online) The phase shifts 4>Ti (t) are shown. The 
TDCIS calculations were done on the one side (red solid) with 
the full electron-ion interaction and all 3d -1 , 4s _1 , 4p _1 ionic 
states active, and on the other side (green dashed) with a 
simplified two-level TDCIS model without residual Coulomb 
interaction (Hi — 0). 



flucncc on the remaining ion via Hi — this means we set 
^residual = q j n p[g_ o, ; wc compare 4>t x (t) obtained from 

the full TDCIS model including Hi and all ionic states 
of Kr + (red solid) and that obtained from the simplified 
two-channel TDCIS model (green dashed) just described. 
Whether or not Hi is included makes no significant differ- 
ence in the behavior of tpT^ (j). Hence, residual Coulomb 
interaction between the ionized electron and the ion has 
almost no effect, i.e., <f>™ sidual (t) w 0. 

3. Field- Induced Mixing with the Neutral Ground State 

The discussions in Sees. HID 1 and HID 2 have shown 
that 0!^ n < 0t and <^ csidual < T . Hence, wc must 
conclude that the main reason for the phase shift comes 
from the field-driven mixing with neutral ground state, 
i.e., 4> t (t) w </4 round (r) [cf. Eq. (22)]. Remember that 
the mixing to the ground state is captured by the terms 
($ | z |$ 4 Q ) , ($?| z |$ ) [cf. Eq. (3)]. These terms are also 
responsible for describing tunnel ionization. To verify 
that the field-induced mixing of the excited TV-electron 
states <I>™ with $o is indeed the main reason for the ob- 
served phase shift, we perform calculations where we once 
switch off the field-driven mixing to 4> after the probe 
step and once where we leave it on. The probe pulse 
is here delta-like, i.e., -E pro be(i) oc S(t — r). We again 
reduce krypton to a two-level atom (excluding the 4s 
orbital) as described in Sec. Ill D 2. In this two-level sys- 
tem, <^ n (r) = <^ osidual (r) = such that the only phase 
shift that can occur is 0|! ound (r) . 

In Fig. 10, the phase 4>t(t) is shown with and with- 



out the mixing to the neutral ground state. If we set 
($o| z |$?) = (i.e., <j% ound = 0) after the probe step, the 
phase shift totally disappears for all pump-probe delays. 
When including the ground state mixing, we obtain again 
the usual behavior of 4>t(t). Hence, we conclude that the 
main source of 4>t(t), which deforms the transition lines 
in the transient absorption spectrum (cf. Fig. 3), is the 
field-induced dressing of the entire iV-electron system — 
particularly the mixing with the neutral ground state. 




-0.5 o 0.5 1 



pump-probe delay % [fs] 

FIG. 10: (color online) The phase shifts 4>t(j) are shown 
for the transition T = ^p^] -1 -»■ pdg^] -1 . The TDCIS 
calculations were done once with (red solid) and once without 
(green dashed) field-driving mixing with the neutral ground 
state <J>o after the probe pulse. In both cases, no dressing of 
the ionic states can occur (see text for details). 
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IV. CONCLUSION 

We have described our theoretical model, namely 
a time-dependent configuration-interaction singles (TD- 
CIS) approach, which has been implemented in the soft- 
ware package xcid. We extended our method and in- 
cluded spin-orbit splitting for the occupied orbitals. This 
extension leads to new symmetry classes and, there- 
fore, new matrix elements. With the help of this multi- 
electron approach, we are able to study attosecond tran- 
sient absorption experiments for overlapping pump and 
probe pulses from first principles. The pump pulse as 
well as the probe pulse are treated non-perturbatively. 

Transient absorption spectroscopy with overlapping 
pulses makes it possible to study the tunnel ionization 
dynamics with a sub-cycle resolution. We find that the 
hole populations extracted from the transient absorption 
spectrum are in close connection to the instantaneous 
hole populations even though this relations is strictly 
speaking only true for non-overlapping pump and probe 
pulses. 

The strong deformations in the transient absorption 
lines, which appear during the ionization process, can 
be fully understood by phase shifts in the induced ionic 
dipole oscillations. We find that the phase shift due to 
the dressing of the ion by the pump pulse is not suffi- 
cient to account for the observed line deformations. We 
also excluded the residual Coulomb interaction between 
the ionized electron and the remaining ion as a possible 
source. 

The main contribution to the phase shift comes from 
field-induced mixing of the excited TV-electron states 
with the neutral ground state, which is highly non- 
perturbative. This dressing mechanism creates a cou- 
pling between the ionized electron and the ion, which 
makes the ionic subsystem dependent on the state of the 



ionized electron and vice versa. This dressing mechanism 
is quite peculiar, since the ionized electron was believed 
to be a spectator, since the probe pulse does only affect 
the ionic states. We are not aware that this effect has 
been discussed in the literature before. 

The non-perturbative mixing to the neutral ground 
state affects also the phase relations between ionic states. 
These phase relations are particularly important for the 
hole dynamics in an atom or molecule [64]. By varying 
the pump-probe delay, the phases between ionic states 
can be influenced, and transient absorption spectroscopy 
provides a way to "read out" these phases. 
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Appendix A: Dipole Polarizability 

The dipole polarizability of a system that is in the 
state S, can be obtained within perturbation theory and 
reads [65] 

„ (q) nSr (S\x m \I) (I\x n \S) 

a m ,„(6)=2> , (Al) 

hii — £js 

where Eg and E] are the energies of the states S and 
J, respectively, and j stands for the sum over all in- 
termediate states. If S is an eigenstate with magnetic 
quantum number Mj, the polarizability entries a x . x (S) 
and a VtV (S) are equal (a X)X (S) = a yty (S)) [66]. The for- 
mulas for ol x x (5) and a z , z (S) read [66]: 



a x . x (S) = A j [J{J-l) + Mj] 



a z , z (S) = 2Aj [J 2 



Mi 



+ Bj [J(J + 1) 
+ IB j Mj 



Mj] 



Cj [(J + 1)(J 
2Cj [(J 



2) 



l) 2 - M 2 



(A2a) 
(A2b) 



where the constants A j,Bj, and Cj stand for specific 
sums over intermediate states / with Jj = Jg — 1, Ji = 
Jg, and Jj = Jg + 1, respectively. 

For singly ionized, atomic krypton, we are interested 
in the polarizabilities for ionic states with a hole in 
the 4p or the 3d orbital manifolds. Using DALTON, 
we calculated polarizabilities that do not include spin- 
orbit coupling. In order to obtain polarizabilities for the 
spin-orbit-coupled ionic states [4p™ J ] _1 and [3c?™ J ] _1 we 
perform first order perturbation theory for degenerate 
states. The diagonal polarizability entries of the spin- 
orbit-coupled ionic states [NLj 1 - 1 ] -1 expressed with the 
polarizabilities of the non-relativistic states \NLml\~ 1 



read 

a n , n (iNLfT 1 ) = £ \ C iZ;m,°? ^ 

a. Ml 

X a n . n ([NLmj.]' 1 ) ■ 

With the help of Eqs. (A2) and (A3) only two polar- 
izability entries a n , n (S),n <G {x,y, z} from states with 
the same N and L quantum numbers are needed and all 
other diagonal polarizability entries of any state S' can 
be obtained as long as S' has the same N and L quan- 
tum numbers. In Table III the calculated polarizabilities 
are shown that are used to obtain all polarizabilities in 
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TABLE III: Static dipole polarizabilities a x>x and a z , z are 
shown for two Kr + states. All other polarizabilities can be 
deduced from these values. Polarizabilities obtained by the 
CIS and CASSCF methods are compared. All values are given 
in atomic units with a precision up to the second digit. 





CIS 




CASSCF 






a Z}Z 




[4P0]" 1 


0.01 


3.13 


10.77 10.46 


[Sdo]- 1 


0.00 


0.00 


9.66 9.52 



Table I. For the CASSCF and CIS calculations we ob- 
tained converged results with the augmented correlation- 
consistent quintuple-zeta basis set including polarization 
functions (aug-cc-pV5Z) [67]. 
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